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ABSTRACT 

Extensive measurements of the X-ray background (XRB) yield a reasonably reli- 
able characterisation of its basic properties. Having resolved most of the cosmic XRB 
into discrete sources, the levels and spectral shapes of its main components can be 
used to probe both the source populations and also alternative cosmological and large- 
scale structure models. Recent observations of clusters seem to provide evidence that 
clusters formed earlier and are more abundant than predicted in the standard ACDM 
model. This motivates interest in alternative models that predict enhanced power on 
cluster scales. We calculate predicted levels and spectra of the superposed emission 
from groups and clusters of galaxies in ACDM and in two viable alternative non- 
Gaussian (x 2 ) and early dark energy models. The predicted levels of the contribution 
of clusters to the XRB in the non-Gaussian models exceed the measured level at low 
energies and levels of the residual XRB in the 2 — 8 keV band; these particular models 
are essentially ruled out. Our work demonstrates the diagnostic value of the integrated 
X-ray emission from clusters, by considering also its dependences on different metal- 
licities, gas and temperature profiles, Galactic absorption, merger scenarios, and on a 
non-thermal pressure component. We also show that the XRB can be used for a upper 
limit for the concentration parameter value. 
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1 INTRODUCTION 

The origin of the cosmic X-ray background (CXRB) is of 
considerable interest due to possible ramifications on the 
properties of the various classes of sources that contribute 
to it, such as clusters and other extragalactic sources. By 
estimating their CXRB contributions and accounting for 
Galactic absorption we can compare with measurements of 
the observed X-ray background (XRB). Since the discovery 
of the XRB (Giacconi et al. 1962) significant effort has been 
devoted to spectrally mapping the emission at energies in 
the range 1CT 1 - 10 2 keV. The radiation has a power-law 
form with a spectral index of T = 1.4 above 1 keV; the spec- 
trum E ■ I{E), where 1(E) is the spectral intensity, peaks 
at E ~ 30 keV and decreases exponentially at higher ener- 
gies. Substantial uncertainties in the measurement of 1(E) 
are due to cross-calibration uncertainties (e.g., De Luce & 
Molendi 2004) and spatial variations reflecting the large- 
scale structure (Gilli et al. 2003; Yang et al. 2003). 



It has been determined that the XRB is mostly the su- 
perposed emissions from discrete extragalactic sources, other 
than at the lowest energies, < 0.25 keV, where its dominant 
contribution is thermal emission from the local bubble (LB) 
and the Galactic halo. Estimated fractional contributions 
of resolved sources to the XRB are 70-80% in the ROSAT 
0.5 - 2 keV (Hasinger et al. 1998) and HEAO-1 A2 2 - 10 
keV (Mushotsky et al. 2000) bands. Other estimates of the 
fractional contributions of resolved sources are 60-80% (Gi- 
acconi et al. 2001), 70% in the 2-8 keV band (Cowie et al. 
2002). More recently, Worsley et al. (2005) determined XRB 
fractions of 85%, 80%, 60%, 50% for the 0.5-2, 2-10, 6-10, 
8-10 keV bands, respectively. Higher fractions were deduced 
from deep Chandra fields (CDFs), 89.57.^% and 86.9^j% 
of the in the 0.5-2 and 2-8 keV bands, respectively (Bauer 
et al. 2004). These contributions include AGN, which dom- 
inate above 1 keV (Hasinger 2004, Brandt & Hasinger 2005 
and references therein; Gilli, Comastri, & Hasinger 2007), 
star forming galaxies, and Galactic stars, but do not in- 
clude emission from high-mass groups and clusters of galax- 
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ies (which will be simply referred to as clusters); see also 
(Hickox & Markevitch 2006,2007). 

A small contribution at very low energies may come 
from the warm hot intergalactic medium (WHIM). Most of 
the unresolved 1-2 keV emission can be accounted for by 
galaxies that are detected in deep HST observations, but 
are too faint to be detected as individual X-ray sources, as 
shown by a stacking analysis (Worsley et al. 2006). Analy- 
sis of the CDFs implies that the maximum contribution of 
clusters is estimated to be ~ 17 % and ~ 20 % in the 0.5-2 
keV and 2-8 keV bands, respectively. Clusters are important 
probes of the primordial density fluctuation spectrum and 
the evolution of the large scale structure. Since conditions 
in clusters can be gauged by levels of X-ray emission, their 
estimated contribution to the CXRB can be used £LS Et St Et _ 
tistical diagnostic measure of the cluster population and its 
evolution. 

The contribution of clusters to the XRB was estimated 
in several works (e.g., Rosati et al. 1998; Gilli, Risaliti,& 
Salvati 1999; Wu & Xue 2001) to be less than ~ 10% of the 
XRB spectral intensity at ~ 1 keV, and even lower at en- 
ergies > 2 keV. The integrated spectral intensity of clusters 
can now be calculated more accurately from our improved 
knowledge of intracluster (IC) gas properties, the relevant 
scaling relations, as well as more precisely determined values 
of the global cosmological parameters. Moreover, X-ray line 
emission, which contributes significantly to the spectrum at 
low energies, particularly in high-mass groups, should be 
added to the thermal Bremsstrahlung emission (as has al- 
ready been shown by Phillips, Ostriker, & Cen 2001). 

Additional motivation for revisiting the issue of clus- 
ter contribution to the CXRB is provided by the possibility 
of using it to constrain alternative cosmological and LSS 
(Large scale structure ) models, that have been proposed 
to remove apparent tension between the predictions of the 
standard model and recent observations of clusters. Possibly 
discrepant results include the following: (a) The detection 
of proto-spheroidal galaxies at z ~ 2 with a substantially 
higher comoving density than what is predicted by most 
semi-analytic models of galaxy formation (Magliocchetti et 
al. 2007). Even more problematic is the detection of struc- 
tures with high velocity dispersions at redshifts z = 4.1 (Mi- 
ley et al. 2004) and z = 2.1 (Kurk et al. 2004). (b) The 'ex- 
cess' CMB power at high multipoles (I w 3000) measured by 
the CBI (Readhead et al. 2004), ACBAR (Kuo et al. 2004), 
and BIMA (Dawson et al. 2006) arrays. The measured power 
levels seem to be too high to be comfortably consistent with 
the predicted CMB primary anisotropy, and are also out- 
side the range of predicted anisotropy induced by the S-Z 
effect (Sadeh, Rephaeli, & Silk 2006, 2007). (c) High val- 
ues of the concentration parameter, which were deduced for 
A1689 from two independent analyses: a joint lensing and 
X-ray analysis (Lemze et al. 2008a, hereafter L08a), and 
a dynamical analysis (Lemze et al. 2008b, hereafter L08b) 
of this cluster. In addition, high concentrations can explain 
the recently measured large Einstein radii in five high-mass 
clusters (Broadhurst & Barkana 2008), obtained from the 
analysis of high-quality ACS and Subaru lensing measure- 
ments and Chandra data. These values prove to be high even 
after relaxing the recent formation time assumption usually 
adopted in theoretical estimations, by considering the prob- 



ability distribution function of halo formation times (Sadeh 
& Rephaeli 2008). 

All the above apparently discrepant observational re- 
sults can be naturally explained in the context of LSS mod- 
els with primordial density fluctuation spectrum that had 
excess power on cluster scales, resulting in earlier formation 
and higher abundances of high-mass clusters than in the 
standard ACDM model. Several such viable models with the 
desired enhanced power on cluster scales are non-Gaussian 
X 2 models (e.g., Koyama, Soda & Taruya 1999) and early 
dark energy models with time-varying coefficient of the 
quintessence equation of state (Doran et al. 2001; Caldwell et 
al. 2003). It has already been suggested by Mathis, Diego, 
& Silk (2004) that scale-dependent non-Gaussianity (at a 
level consistent with observational constraints from WMAP 
measurements) may be probed through its detectable mani- 
festations on clusters. The predicted S-Z power spectra and 
cluster number counts were determined by Sadeh, Rephaeli, 
& Silk (2006, 2007). 

Earlier formation, higher abundances, and more con- 
centrated IC gas profiles will also be reflected in possibly 
discernible differences in their X-ray properties when com- 
pared with those predicted in the standard ACDM model. 
Therefore, in addition to the need to update and improve 
calculations of the contribution of clusters to the CXRB, 
further motivation for extending our previous work on alter- 
native cosmological models is provided by the expectation 
that cluster X-ray properties could yield additional useful 
constraints on these models. 

This paper is arranged as follows: In §[2] we describe the 
methodology adopted for the calculation of the spectral con- 
tribution of clusters and high mass groups in three distinct 
cosmological models. We present our results in § [3l where 
predicted XRB levels are made for alternative cosmological 
models, different levels of Galactic absorption, various gas 
profiles, cluster mergers, and in the case when the total gas 
pressure includes a non-thermal component. We conclude 
with a discussion in § 0] 



2 METHODOLOGY 

X-ray emission from clusters is calculated in ACDM, non- 
Gaussian (NG), and early dark energy (EDE) models. In the 
NG models the probability distribution function (PDF) of 
the primordial density fluctuation field has a Xm form. The 
largest deviation from Gaussianity occurs when there is one 
such field (m = 1). We consider the m = 1 and m — 2 cases. 
A full account of the the non-Gaussian and early dark energy 
models that we have used are provided in Sadeh, Rephaeli, 
& Silk (2006, 2007). Here we present a brief summary of the 
most important aspects which are of particular relevance to 
this work. 



2.1 ACDM and Non-Gaussian Models 

The number density of clusters as a function of mass and 
redshift is calculated using the Press-Schechter (Press & 
Schechter 1974, hereafter PS) mass function 

dn(M) p m da 
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The top and bottom expressions refer to the Gaussian and 
models, respectively, and TV is a normalization factor, 
introduced in order to arrange that all the matter in the 
universe will be incorporated into halos. 

For the Gaussian, NG (m=l), and NG (m=2) mod- 
els, N = 2,0.31,0.37, respectively. Additional quantities 
appearing in the PS mass function include v = ^pjy, 
which denotes the critical density for collapse in units 
of the mass variance. The latter quantity is defined as 
a 2 (M) = j^jj J d 3 kW% H (kR)P(k), where W TH (kR) = 
(fcllp [sin(kR) — kRcos(kR)] is the Fourier transform of the 
Top-Hat window function, and the relation between M and 

R is R = ^ 4^- ) • The power spectrum is taken as 

P(k) = A p k n T 2 (k), with n = 1 and n = -1.8 for the Gaus- 
sian and NG models, respectively. 

Adiabatic CDM and isocurvature transfer functions are 
used in the Gaussian and NG models: 

M 2 + 34 g 34g? I 1 + 3 - S9q + C 16 - 1 ?) 3 + (5-46g) 3 
+ (6.71g) 4 ]- ' 25 
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(3) 



where q = fe/(fi m h 2 )Mpc 1 (Bardeen et al. 1986), and A p 
is the normalization set by erg: 



(4) 



The cosmological parameters were taken from the WMAP 
3-year data, with Qa = 0.76, fl m = 0.24, h = 0.73, and 
erg = 0.74 (Spergel et al. 2007). Due to the steep dependence 
of the cluster mass function on erg, we have accounted for 
the higher value, ~ 0.82, inferred from WMAP5 (Komatsu 
et al. 2008). This point is elaborated on in § [4] The gas 
profile is taken to be a (5 model, i.e. n ga s(r) = n Bas (0)(l + 
( — ) 2 ) -1,5 ' 3 , where the central density is normalized such 
that the total gas mass constitutes a fraction 0.1 of the total 
cluster mass included within the virial radius, and r c denotes 
the core radius. To allow for a temperature structure, the 
profile is taken to have a polytopic form, T(r) = T(0)(1 + 
(~ ) 2 ) -1 ' s ' 3 (" T— l ', where 7 is the polytopic index, and the 
central temperature can be derived under the assumption of 
hydrostatic equilibrium (HE), that yields 



T(r = 0) = 



\},m v G M v i r 
SP'ykB r vir 



1.5/3( 7 -l) + l 



(5) 

where G, m p , and ks are the gravitational constant, pro- 
ton mass, mean molecular weight, and the Boltzmann con- 
stant, respectively. 

We first consider isothermal gas, 7 = 1, and take 



Ca 



10, 



|. Note that C gaa is the virial ra- 



dius in units of the gas core radius, essentially the analogous 
concentration parameter, where the virial radius is in units 



of the NFW (Navarro, Frenk, & White 1997) scale radius. 
We consider (in § 13 -3 P a range of values of C go ,s and 7 in 
order to assess the dependence on these parameters. With 
the gas density and temperature profiles specified, we can 
calculate the spectral luminosity of a cluster with mass M 
and redshift 2 



-(E) = 



EECtV, 



(6) 



where ee is the specific (spectral) emissivity in 
erg sec -1 cm~ 3 kev~ x . We compute the emissivity 
with the MEKAL plasma code and compare with results 
obtained when the emissivity includes only Bremsstrahlung 
continuum (Fig. [TJ. The calculation of the Gaunt factor 
was carried out using the analytic fitting formula of Itoh et 
al. (2000). 

We consider all systems in the mass range 1.1 • 10 13 — 
7.3 ■ 10 15 h^ M Q , which roughly includes all halos larger 
than a cD galaxy. Increasing the upper halo mass limit does 
not affect the results, since the PS formalism predicts an 
extremely low abundance of such massive halos. Indeed, the 
most massive currently known cluster, A370, has a mass 
(2.93±o:! • !0 15 K.t M ©> Broadhurst et al. 2008) which is 
less than half the assumed upper end of this range. The 
specific flux from all clusters in this mass range is 



Jo JM min dMdz 4nD 2 L (z) 

(7) 

where Dl is the luminosity distance; in our calculations we 
take z max — 3. 



1 



2.2 EDE Model 

The evolutionary history in the EDE model is different from 
that in the ACDM universe, for which we used the Gaussian 
and NG models. The vacuum (Quintessence) energy density, 
Pq(z), is redshift dependent, 



PQ = PQo(l + *) 3l ™ ] , 



(8) 



where uj(z) is the equation of state coefficient ( — 1 for a 
cosmological constant). In our calculations we adopt the 
parametrisation of Wetterich (2004): 



Uj(z) 



l+ulog(l + z) 



(9) 



and 



log - 



i+log 



where fl e is the early 



quintessence density and ujo is uo{z = 0). The specific EDE 
model assumed here is that investigated by Bartelmann, Do- 
ran, & Wetterich (2006; hearafter BDW), for which n e = 
8 • 10" 4 and cu = -0.99. 

Several cosmological quantities which are affected by 
the presence of early dark energy require re-evaluation. 
These include the (comoving) radial and angular diameter 
distances, linear growth factor, the critical density for spher- 
ical collapse, and the overdensity at virialization, all of which 
are explicitly calculated by Sadeh, Rephaeli, & Silk (2007). 
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Figure 1. The integrated contribution of clusters and high mass 
groups with masses (1.1 ■ 10 13 - 7.3 ■ 10 15 ) hp \ Mq to the 
XRB intensity as a function of energy is shown together with 
the XRB data. Emission processes include free-free (black dash- 
dotted curve), optically thin hot plasma with zero metal abun- 
dance, A=0 (red dashed curve), and optically thin hot plasma 
with A=0.3 (blue solid curve). 



3 RESULTS 

3.1 Cluster contribution to the CXRB in ACDM 

As noted, we have used the MEKAL code to calculate 
gas emissivity. Inclusion of only Bremsstrahlung contin- 
uum emission appreciably underestimates radiation from 
hot plasmas, especially at the relatively low temperatures of 
the more abundant groups and poor clusters, for which line 
emission at low energies and continuum free-bound emission 
processes should also be included. In Fig. [T] we show the su- 
perposed spectrum of clusters in the mass range 1.1 • 10 13 
\ Mq - 7.3 • 10 15 Mq when including only free-free 
emissivity (dash-dotted line) , optically thin hot plasma with 
zero metal abundances, A=0 (dashed line), and optically 
thin hot plasma with A=0.3 (solid line). The XRB measure- 
ments (with uncertainties), shown here by the green region, 
are based on the compilation of Gilli et al. (2007) , who com- 
bined measurements from HEAO-1 (Gruber 1992; Gruber et 
al. 1999), ASCA GIS (Kushino et al. 2002), ROSAT PSPC 
(Georgantopoulos et al. 1996), two different sets of XMM 
measurements (Lumb et al. 2002; De Luca & Molendi 2004), 
ASCA SIS (Gendreau et al. 1995), BeppoSAX (Vecchi et al. 
1999), RXTE (Revnivtsev et al. 2003), and shadowing ex- 
periments (Warwick & Roberts 1998). These measurements 
differ slightly in their minimum Galactic latitude. The effect 
of different minimum Galactic latitudes on the spectrum is 
investigated in § 13.21 As can be inferred from Fig.fj] the level 
of the Bremsstrahlung emission alone is quite low in compar- 
ison with the XRB measurements. Its contribution peaks at 
~ 0.4 keV, where it constitutes ~ 20% of the CXRB. Inclu- 
sion of the full continuum (Bremsstrahlung and free-bound) 
with A=0 enhances the cluster emissivity by ~ 1.25 across 
the entire 0.1 — 10 keV energy band. With the higher abun- 
dance of A=0.3, the emissivity is much higher, especially in 
the 0.4 — 1 keV energy band. 

Clearly, high mass groups and clusters contribute to 
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Figure 2. The relative XRB spectral intensity contributed by 
high-mass groups and clusters, defined by the respective mass 
ranges of 1.1 ■ 10 13 h~* M - 10 14 h~J, Mq and 10 14 h~\ Mq - 
7.3 • 10 15 h~J, Mq. 

the CXRB in different energy bands. Clusters contribute 
at higher energies due to their higher masses. We illustrate 
in Fig. [2] the contribution of the two populations, where 
we took the mass range of the high mass groups and clus- 
ters to be 1.1 • 10 13 hu"? M - 10 14 h^i, Mq and 10 14 h^ 
Mq - 7.3 • 10 15 Mq, respectively. The former popula- 
tion dominates the emission below ~ 2.5 keV, whereas the 
latter contributes mostly at higher energies. Clearly, the su- 
perposed emission of the lower mass systems can be used to 
constrain models that predict large spectral bumps, whereas 
the higher mass clusters can be used to constrain models in 
which fractional contribution of clusters is higher than the 
low residual level, namely upon subtraction of the dominant 
AGN contribution. 

We have explored the dependence of the predicted in- 
tensity on the various LSS quantities and IC gas model. Here 
we consider the impact of using approximate fitting formu- 
lae for the linear growth factor, critical density for spher- 
ical collapse, and overdensity at virialization. To do so we 
numerically solved the relevant equations describing these 
quantities (see, e.g., Sadeh, Rephaeli, & Silk 2007) and used 
these exact solutions instead of the respective approximate 
fitting formulae. We find that at low energies the two sets 
of calculations yield essentially the same results. However, 
with increasing energy the exact calculation yields a level of 
intensity which is lower by as much as ~ 17% at 10fce^,as 
can be seen in Fig. 

In the following subsections we assess the impact of IC 
gas modeling uncertainties, merger histories of clusters, and 
dependence of the predicted CXRB intensity on Galactic 
absorption. 

3.2 Galactic absorption 

To account for photoelectric absorption in the Galaxy, we 
combined our calculated CXRB intensity with the XSPEC 
WABS model. We then computed the average column den- 
sity in the following way: 10 4 lines of sight (los) were picked 
at random with a uniform probability distribution on a 
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Figure 3. The spectral intensity calculated using analytical fits 
for the linear growth factor, critical density for spherical collapse, 
and overdensity at virialization, compared with the corresponding 
spectral intensity calculated using the exact numerical solutions 
to the differential equations for these quantities. (This comparison 
was made for the case A = 0.) 

sphere. This number proved to be sufficient in order to pro- 
duce statistically meaningful results by virtue of the fact 
that with 10 3 los the average column density changed negli- 
gibly and its error was still small. Note that the number of 
los has an upper limit of 41259 since the resolution of the 
column density is ~ 1°. 

For each los we obtained the column density using the 
nh tool (Dickey & Lockman 1990). This procedure was re- 
peated 10 times, yielding the average column density and 
its standard deviation, which was taken to be the error. 
In this way we calculated the average column density of 
the Galaxy in all directions. In order to take into account 
the fact that XRB observations are performed at different 
Galactic latitudes, we calculated the average Galactic ab- 
sorption at stripes of Galactic latitude b > 10° and b > 30°, 
namely excluding the disk region, and the combined regions 
of the disk and the bulge, respectively. In addition, many 
X-ray observations point at regions on the sky where the 
column density is significantly lower, the goal of which is 
the detection of the soft X-ray background (which is ab- 
sorbed in regions of high column densities). Such a region 
is the Lockman Hole, associated with a minimum hydrogen 
column density Nh = (4.5 ± 0.5) • 10 19 cm -2 (Lockman et 
al. 1986). 

Results for the average column density at different re- 
gions marked by the lower limit of the Galactic latitude are 
summarized in table [1] The contribution of clusters to the 
XRB with A = 0.3, taking into account the effect of Galactic 
absorption, but excluding the effect of mergers, is illustrated 
in Fig. |4] 

3.3 IC Gas Density and Temperature Profiles 

The density profile of IC gas is usually represented by a (3 
model, which fits well X-ray and SZ measurements (though 
occasionally a double /3 model provides a better fit). How- 
ever, the assumption of isothermality is known to be unreal- 



Table 1. The average column density at latitudes higher than 
b. The Lockman hole is situated at a = 10 h 45 m , <5 = 57°2o' 
(Lockman et al. 1986) 



b n 


(N H (> 6)> [10 20 cm" 2 ] 





13.9 ±0.3 
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6.05 ± 0.06 
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3.54 ±0.04 
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Figure 4. The impact of Galactic absorption on the integrated 
contribution of clusters and high mass groups to the XRB spectral 
intensity, for Galactic latitudes b > 0° (purple striped curve), 
b > 10° (black dash-dotted curve), b > 30° (red dashed curve), 
and the Lockman hole (blue solid curve) 

istic, especially in the outer regions of clusters. A polytropic 
equation of state provides a more realistic temperature pro- 
file, since at r » r c , n(r) oc r~ 313 and T(r) oc r ~ 3,3(7 " 1) . 

We can use insight on the profile of the gas entropy 
parameter, K = Tn~ 2 ^ , to select the relevant range of 7 
values. In the outer regions of clusters, r > O.liivir, beyond 
the cool core, the entropy is entirely due to gravitational 
processes, and its profile is expected to be a featureless power 
law approaching K oc r 1 ' 1 (Tozzi & Norman 2001; Voit, 
Kay, & Bryan 2005). However, a recent model-independent 
analysis of X-ray and lensing measurements of A1689 yielded 
0.82 ± 0.02 for the power- law index of the profile (L08a). 
For isothermal gas, /3 = 2/3 implies K oc r 4 ^ 3 . This value 
of the index is higher than the two aforementioned values 
of 1.1 and 0.82. Thus, in addition to not being realistic at 
large radii, the isothermal case is also inconsistent with the 
entropy profile when j3 — 2/3 is adopted. In order for the 
entropy profile to have an index of 0.82, 7 has to be ~ 1.25, 
which is in good agreement with both ASCA measurements 
of 30 nearby clusters (Markevitch et al. 1998) and numerical 
simulations (Loken et al. 2002). 

In fig. [S] we explore the cluster contribution to the 
CXRB for different values of C ga s and 7. For C gas = 10 
and 7 = 1.25 the X-ray intensity is higher at E > 1 keV as 
compared with its value for 7 = 1, reflecting the higher cen- 
tral temperature. Although the temperature at r > r V i r is 
lower than in the 7 = 1 case, the radial integration (to r V i r ) 
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Figure 5. Dependence of the cluster CXRB contribution on dif- 
ferent values of C ga s and 7. 



only results in a boost of the intensity level at high energies. 
For C gas = 5, there is no significant difference between the 
7=1 and 7 = 1.25 cases, a result of which can be explained 
taking into account the fact that the change in the central 
temperature (which depends on C gaa ) is lower, and that the 
core radius of the gas is bigger so that the range of difference 
between the core and the virial radii is smaller. 



3.4 Impact of cluster mergers 

The formation of structure in CDM models is hierarchical; 
clusters evolve through mergers of subclumps. During these 
merger events the temperature and luminosity of IC gas are 
temporarily boosted up with respect to the corresponding 
equilibrium values, in a manner that has been characterized 
by hydrodynamical simulations. Neglecting this boost may 
lead to an underestimation of the cluster contribution to the 
CXRB. Randall, Sarazin, & Ricker (2002, hereafter RSR02) 
employed a semianalytic model to estimate the effect of 
merger boosts on the X-ray luminosity function (XLF) and 
temperature function (TF). A quantification of the luminos- 
ity and temperature spikes in individual mergers was made 
based on N-body /hydrodynamical simulations, whereas the 
statistics of the merger history were determined from Ex- 
tended PS (EPS) merger trees. RSR02 used these merger- 
boosted TFs and XLFs to infer as and fi m by converting 
them to PS mass functions and comparing with local (z — 0) 
and distant (either z = 0.5 or z = 1) cluster samples. For 
the ACDM model an increment of 20% in erg was deduced by 
fitting to both the TF and XLF. On the other hand, results 
for fi m were inconclusive; fitting to the TF and XLF ex- 
hibited a decrement and an increment, respectively. In fact, 
the results obtained from the XLF fit were less reliable, as 
indicated by the fact that fi m inferred from the non-boosted 
functions was lower by 60% than the actual value used to 
construct the merger trees. The reason for the lower degree 
of reliability associated with the fit to the XLF is due to the 
fact that the range of X-ray luminosities is broader than that 
of X-ray temperatures; put otherwise, the X-ray luminosity- 
temperature scaling is steeper than a linear relation. As a 
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Figure 6. The cluster contribution to the CXRB including the 
effect of mergers. Plotted are the XRB data (green area), the clus- 
ter contribution without the effect of mergers (blue solid curve), 
and the corresponding contributions when their effect was mod- 
elled by selecting values of as higher by 10% and 20% (red dashed 
and black dot-dashed curves, respectively), and by increasing as 
by 20% but lowering Q m by 20% (purple stripped curve). Also 
shown is the currently estimated upper limit on the residual (i.e., 
after accounting for emission from AGN) 2-8 keV emission and 
its 1-a range (cyan solid and dashed lines, respectively). 

result, the overall effect of luminosity boosts on the XLF is 
not as strong as the effect of temperature boosts on the TF. 

The fitted er 8 to the boosted TFs was higher by 20% 
than the value used to construct the trees, whereas with the 
non-boosted TFs it was higher by only 10% (a difference 
which is due to the approximation used to determine the 
cluster formation redshift). In Fig. [6] we compare our calcu- 
lated cluster contribution to the CXRB including the effect 
of mergers with the data. The impact of mergers was mod- 
eled in three different ways by selecting values of as higher 
by 10% (as = 0.81) and 20% (rr 8 = 0.89), and by adopting 
a combination of as higher by 20% (as = 0.89) and of Q m 
lower by 20% (Q m — 0.19). As is clear from the figure, the 
third choice, i.e. the one combining a higher as and a lower 
Q m , would seem more likely, simply by virtue of its consis- 
tency with the data. In fact, this scenario was also preferred 
by RSR02, who noted that the effect of merger boosts on 
as and Q m depends somewhat on the detailed method used 
to determine the TFs and XLFs, and the criteria used to fit 
them. It should be noted that this holds for C ga s ~ 10; for 
C 9 as = 5, none of the models exceed the observed intensity 
range. However, this is no longer the case if as is higher, as 
is alluded to in § U In addition, in the the 2 — 8 keV energy 
band high values of as boost the cluster CXRB contribution 
above the 1-er upper limit of the residual XRB, which was 
calculated by multiplying the measured XRB upper bound 
(in the 2 — 8 keV band) by the residual fraction determined 
by Bauer et al. (2004). 

Results for merger boosts that emerge from various 
works differ in several respects. Bullock et al. (2001) found in 
their simulations a value of \ m ed, the spin parameter (Pee- 
bles 1971), lower than 0.05, the value used in RSR02. This 
would slightly increase the effect of merger boosts. On the 
other hand, preheating and non-gravitational effects may 
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Figure 7. The cluster CXRB contribution for a purely thermal, 
and for the case when the thermal pressure is only 90 %, or 80% 
(blue solid, red dashed, black dot-dashes curves, respectively). 



somewhat reduce the effect of merger boosts since the frac- 
tional increase in the X-ray luminosity and temperature 
would be smaller if non-gravitational heating is considered. 
This effect would dominate in low-mass clusters (Ponman et 
al. 1996). 



3.5 Non-thermal pressure 

The hydrostatic equation (HE) for an ideal gas can be inte- 
grated from a given radius r (assuming spherical symmetry) 
out to some limiting radius 

GM(^ r')fj,m p pgas(r')dr' 



[p 9 a 3 (r)T(r)] = C 



(10) 



where p gas and T are the gas mass density and temperature, 
respectively, M(^ r) is the total mass within radius r. If the 
only contribution to the pressure - which balances gravity - 
is the thermal gas, C = 1 (e.g., Binney & Tremaine 1987). 

In our work so far we have assumed that clusters are 
in hydrostatic equilibrium, and that gas pressure is strictly 
thermal. However, this may not necessarily be the case, at 
least in some clusters. Accordingly, the impact of having a 
second pressure component was assessed heuristically (ig- 
noring other possible ramifications of the source of the ad- 
ditional pressure), simply by gauging the lower IC gas con- 
tribution by taking C = 0.9 and C = 0.8 (assuming, rather 
realistically, that a non-thermal - such as turbulent and en- 
ergetic particle - component does not exceed 20% of the 
total pressure). We can see from Fig. [7] that as C decreases 
the cluster contribution is lowered at energies > 1 keV, with 
practically no change at lower energies. 

3.6 Seth & Tormen mass function 

The PS mass function, which is derived under the naive 
assumption of spherical collapse, underpredicts the abun- 
dance of higher mass clusters, and overpredicts the abun- 
dance of lower mass systems, with respect to those observed 
in N-body simulations. Such simulations manifest a clear 
non-spherical collapse behaviour. Sheth & Tormen (1999, 



Figure 8. The cluster CXRB contribution calculated with the PS 
(solid line), ST (dashed line) mass functions; also shown is the re- 
sult for ST mass function with C ga3 = 5 and Galactic absorption 
with b m i n = 10° (dash-dotted line). Also shown is the currently 
estimated upper limit on the residual (i.e., after accounting for 
emission from AGN) 2-8 keV emission and its 1-cr range (cyan 
solid and dashed lines, respectively). 



hereafter ST; see also Maggiore & Riotto 2009) developed a 
modified mass function, which provides a better agreement 
with results of simulations (Sheth, Mo, & Tormen 2001): 



F(v) 



1 + 



,,2 \ v 
a 



(11) 



where A = 0.3222, a = 0.707, and p = 0.3. The PS mass 
function is reproduced when A — 0.5, a — 1, and p = 0. 
In Fig. [8] we plot the CXRB spectrum in the mass range 
1.1 • 10 13 ho.7 M - 7.3 • 10 15 M , using the PS (solid 
line) and ST (dashed line) mass functions, and the ST mass 
function with Cgas ~ 5 and Galactic absorption with b m ; n = 
10° (dash-dotted line). 

Since the ST mass function generates a higher popula- 
tion of high-mass halos with respect to the PS mass function 
in almost the entire considered mass range, it gives rise to a 
significant increase of the predicted spectrum. If we neglect 
Galactic absorption, the computed levels exceed the mea- 
sured XRB range; if Galactic absorption is included, the 
power spectrum peaks at approximately the same level as 
the unabsorbed spectrum generated with the PS mass func- 
tion. Since not all of the XRB in the 0.4 < E < 1 keV energy 
range is due to clusters, their partial contribution actually 
puts a very strong constraint on the value of Cgas'- As is clear 
from the figure, the ST mass function yields intensity lev- 
els which are roughly half of the observational upper limit, 
if Galactic absorption with b m i n — 10° and C ga s = 5 are 
adopted. Note that the residual emission in the 2 — 8 keV 
band yields an upper limit of ~ 10 on C ga s- 



3.7 Alternative cosmological models 

As discussed in the § [1] there is considerable interest in vi- 
able alternative cosmological models that predict enhanced 
power on the scale of clusters, and their earlier formation 
than in ACDM. In this section we determine the CXRB 
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Figure 9. X-ray spectral intensity curves produced by a cluster 
population derived from a PS mass function with a Gaussian 
(solid line), \\ (dashed line), and x| (dotted line) PDFs, and 
an EDE model (dot-dashed line). Also shown is the currently 
estimated upper limit on the residual (i.e., after accounting for 
emission from AGN) 2-8 keV emission and its 1-rr range (cyan 
solid and dashed lines, respectively). 



cluster contribution in the NG and EDE models specified 
above. In Fig. [9] we compare the X-ray intensity spectra 
produced by a cluster population derived from a PS mass 
function with a Gaussian (solid line), \i (dashed line), and 
xl (dotted line) PDFs, and an EDE model (dot-dashed line). 
As is clear from Fig. [9] the predicted intensity levels in the 
NG models for both m=l and m=2 are much higher than 
the observed range below ~ 1 keV. Clusters overproduce the 
XRB at a level which is well beyond the range of modeling 
uncertainties. The predicted level of contribution is in fact 
excessive also in the 0.5 — 2 keV band, where clusters are pre- 
sumed to contribute at most 10-20% of the measured XRB 
intensity. In the 2 — 8 keV energy band the predicted in- 
tensity levels in these models are higher than the 1-a upper 
limit of the residual XRB. This interesting result strengthens 
the conclusion of Sadeh, Rephaeli & Silk (2007) who showed 
that the level of CMB anisotropy induced by the S-Z effect 
is significantly higher than the current upper limits set by 
the BIMA experiment. Based on these excess intensity and 
power levels, these particular NG models seem to be ruled 
out. 

The degree at which predictions of the EDE model may 
be reconciled with the observational constraints depends on 
the normalization of the matter power spectrum. With the 
same as as derived from WMAP, the intensity bump exceeds 
the XRB data, similar as in the non-Gaussian models. How- 
ever, when we determine the value of as by requiring that the 
cumulative halo number at z = in this model matches the 
corresponding value in the A CDM we obtain essentially the 
same spectrum. The predictions for the surface brightness 
distribution as a function of mass and redshift are shown in 
Fig. QIJ]for the three models. 

The level of tension between predictions of the EDE 
model and the CXRB measurements depends very much 
on the value of as- Normalizing as such that the cumula- 
tive halo number at z = in this model is roughly equal 
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Figure 10. The redshift (top panel) and mass (bottom panel) 
distributions of the CXRB spectral intensity contributed by clus- 
ters are shown for the Gaussian (blue solid curve), x\ ( re d dashed 
curve), x| (purple dotted curve) PDFs, and for the EDE model 
(black dot-dashed curve). 



to the corresponding value in ACDM leads to similar pre- 
dicted spectra. This conclusion is valid even if we adopt 
Cgas = 5, in addition to using the WMAP5-deduced as, 
which is higher by ~ 10 % than the corresponding WMAP3 
value (adopted here). The impact of the higher value of as 
is illustrated in Fig. [6] 



4 DISCUSSION 

Significant advancements in sensitivity and spatial resolu- 
tion have resulted in an improved ability to resolve the XRB 
into its different components and sources. The combination 
of high quality data and enhanced-accuracy calculations can 
be used in order to constrain different aspects which affect 
directly and indirectly the XRB levels. As elaborated upon 
in § 13.11 the spectral bump in the band 0.4 < E < 1 keV, 
an important characteristic of line emission from high-mass 
groups and poor clusters, constitutes a particularly sensi- 
tive probe; models which boost the emission in this energy 
band might exceed the XRB data. Another important ob- 
servational constraint includes the XRB levels in the hard 
energy band, 2 — 10 keV. It seems to have been established 
that AGN constitute the major contributors in this energy 



Cluster contribution to CXRB 9 



band. Consequently, models in which the cluster contribu- 
tion is appreciably higher than ~ 10% — 20% of the total 
CXRB intensity would likely be non-viable. Moreover, the 
combination of extensive measurements of the CXRB and 
detailed predictions of the contributions of AGN, clusters, 
and starburst galaxies, provide a diagnostic tool of the num- 
ber density and evolution of these sources, as well as a probe 
to set constraints on alternative cosmological models. 

We demonstrated in this paper that meaningful con- 
straints can be set on alternative cosmological models by 
contrasting the predicted contribution of clusters to the 
CXRB in non-Gaussian Xm an d EDE models with the mea- 
surements. As discussed in § [TJ there seems to be some ob- 
servational evidence for earlier cluster formation and higher 
abundances (as compared with those in ACDM) that are 
predicted in these models. Our basic result is that the clus- 
ter contribution in the Xi an d X2 models is excessively high 
at low (< 1 keV) energies, so much so that these models 
could probably be ruled out. Moreover, in the high energy 
band, 2 — 10 keV, the predicted spectral intensity attains 
unrealistically high levels by virtue of the fact that most of 
the CXRB at these energies originates in AGN. The viabil- 
ity of the EDE model considered here depends very much 
on how exactly this model is normalized with respect to the 
standard model. With the same as as derived from WMAP5 
values, the predicted energy bump in this EDE model ex- 
ceeds the measured range, as is the case in the non-Gaussian 
models. However, normalizing as such that the cumulative 
halo number at z — in this model matches the correspond- 
ing value in ACDM, yields nearly the same spectrum. These 
conclusions are still valid even if we adopt C ga s = 5 and 
Galactic absorption with b m i n = 10°, along with ST mass 
function and the value of as deduced from WMAP5 mea- 
surements, which is higher by ~ 10 % than the correspond- 
ing WMAP3 value. The effect of this increment in as is 
illustrated in Fig. [S] 

We also showed how the CXRB can be used to con- 
strain merger scenarios. We estimated the impact of merg- 
ers heuristically by determining the resulting boost in the 
CXRB in terms of the change in as and fl m , following the 
work of RSR02. CXRB levels calculated for two of the three 
models exceed the observational data in the spectral bump 
region. The 20% increment of erg leads to an overproduction 
of the XRB intensity levels, also when adopting C ga s = 5 
and the WMAP5 parameters, whereas a 10% boost is in 
marginal agreement with the data (for a more detailed ac- 
count of the models see §[S]). The impact of Galactic absorp- 
tion, neglected in these calculations, is discussed in § 13.21 

We explored the impact of non-gravitational heating on 
the CXRB and found only a minor change. While the spec- 
tral bump is not affected at all, higher energy emission from 
clusters decreases in proportion with the amount of non- 
gravitational heating. Nonetheless, this decrement is small, 
unless non-gravitational heating is considerable. Changing 
the polytropic index, 7, from 1 (isothermal) to 1.25 also 
does not alter the spectrum considerably, especially for low 
C g as- With Cgas = 10 and 7 = 1.25 the central temperature 
is higher and the temperature at large radii, r > r V i r is lower 
than with the 7 = 1 case. However, since the integration of 
the mass and gas profiles is performed up to the virial ra- 
dius, we only notice an enhancement of the spectrum levels 



at high energies. For C ga s = 5 there is almost no difference 
between the 7 = 1 and 7 = 1.25 cases. 

It was shown that inclusion of all relevant emission 
processes (i.e., not just Bremsstrahlung) is very important. 
Especially important is line emission, which generates the 
0.4 < E < 1 keV spectral bump. We calculated the spec- 
trum taking mean abundances of and 0.3 solar. In clusters 
the metallicity is typically ~ 0.2 solar outside the central 
region (Baldi et al. 2007; Leccardi & Molendi 2008). The 
Fe abundance is nearly similar in all hot clusters, > 4.5 
keV (Mushotzky 2004 and reference within), even though 
medium temperature clusters, 2 — 4 keV, have abundances of 
~ 0.4 (Baumgartner et al. 2005). There seems to be a slight 
redshift dependence; XMM and Chandra show no evolution 
in Fe abundances up to z w 0.8 (e.g., Mushotzky 2004 and 
reference within). However, a certain degree of enrichment 
must have occurred. Balestra et al. (2007) found that the 
average iron content of IC gas at z = is a factor ~ 2 
higher than at z =~ 1.2. Several low-redshift, bright X-ray 
groups were shown by Finoguenov & Ponman (1999) to have 
Fe abundances of < 0.3 at radii higher than ~ 100 kpc. For 
almost any kind of abundance gradient and evolution the 
CXRB contributions of high mass groups and clusters is be- 
tween our derived values for A — and A = 0.3. 

The impact of using a mass function different than PS 
was explored by repeating the calculations with the ST mass 
function, which is characterised by higher abundances of 
high- mass halos in the mass range M > 4-10 13 h^"* Mq with 
respect to the corresponding numbers predicted by the PS 
mass function. The enhanced abundances result in apprecia- 
bly higher intensity levels that possibly exceed the measured 
XRB range for high values of C ga s even in ACDM. Clearly, 
given that not all of the XRB in the 0.4 < E < 1 keV 
energy range is due to high-mass groups and clusters, this 
implied excess provides a very tight constraint on the maxi- 
mum value of Cgas ■ (Note that the mass range of the Jenkins 
et al. (2001) mass function, which was deduced from N-body 
simulations, is irrelevant to our work here.) If C ga s ~ Cdm, 
where Cdm is the concentration parameter, the XRB can 
put strong constrains on the value of the concentration pa- 
rameter. 

Finally, we note that whereas we have accounted for 
emission from groups and clusters, our treatment does not 
include emission from the filamentary WHIM. With temper- 
atures ~ 10 or slightly higher, WHIM emission could also 
contribute somewhat to the CXRB. Clearly, any additional 
contribution will only strengthen our conclusions on the vi- 
ability of alternative models in which the predicted spectral 
bump exceeds the observed range. 
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